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Abstract 

We propose and study a new one-dimensional traffic flow cellular automaton 

(CA) model of high speed vehicles with the Fukui-Ishibashi-type acceleration 

for all cars and the Nagel-Schreckenberg-type (NS) stochastic delay only for 

the cars following the trail of the car ahead. The main difference in the 

delay scenario between the new model and the NS model is that a car with 

spacing ahead longer than the velocity limit M may not be delayed in the new 

model. By using a car-oriented mean field theory, we derive the fundamental 

diagrams of the average speed as the function of car density analytically. Our 

theoretical results are in excellent agreement with numerical simulations. 
PACS numbers: 64.60.Ak, 05.40. -Fj, 05.70.Jk, 89.40. -Fk 
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I. INTRODUCTION 



Traffic fiow cellular automatoii(CA) models have attracted much interest recently. Com- 
pared with other dynamical approaches, e.g. the fluid dynamical approach, to this class of 
problems, CA models are conceptually simpler and can be easily implemented on computers 
for numerical investigations [1-4]. 

Two popular one-dimensional (ID) traffic flow models are the Pukui-Ishibashi(FI) model 
[5] and the Nagel-Schreckenberg (NS) model [6]. An exact car-oriented mean- field (COMF) 
theory has been developed for the FI model with arbitrary limit on the maximum speed 
t^max, car density p and delay probability / [7,8]. However, for the NS model with high 
speed vehicles (vmax > 1) and stochastic delay, there is still no established exact analytical 
theory up to now [9,10]. 

The acceleration and stochastic delay rules of the NS model lead to comphcations in the 
time evolution of the flow and hence it is very difficult for exact analytical studies. In order 
to understand how these rules affect the evolution and the corresponding asymptotic state, 
we study a new ID traffic flow CA model in which only the cars following the trail of the 
car ahead may be delayed. 

The plan of the present paper is as follows. The definition of the model and the evolution 
equations for the inter-car spacings are given in Sec. II. In Sec. Ill some observations are 
made to describe the steady state of the system. We present the fundamental diagrams 
for the low car density case with arbitrary vehicle speed limit and the high density case 
with vehicle speeds hmited to 1 and 2 in Sec. IV. Excellent agreements between numerical 
simulations and theoretical results are shown in Sec. V together with a discussion on our 
results in connection to the FI and NS models. 



2 



II. THE MODEL 



Let N be the total number of cars on a ID road of length L. The density of cars is 
p = N/L. Let Cn{t) be the number of empty sites in front of the nth car at time t, and 
Vn{t) be the number of sites that the nth car moves during the time t step. 

The new model adopts the following acceleration rule [5]: 



We call the nth car "a car that follows the trail of the car ahead" if = Cn{t). It 

means that the nth car may become the neighbor of the car ahead if the car in front stops. 
Stochastic delay is introduced in such a way that all the cars which follow the trail of their 
cars ahead have a probability / to move forward one site less than it is allowed by Step 1, 
i.e., we have 

Step 2: Vn{t) = v'^(t) - 1 with the probability /, if v'^{t) = Cn{t) and v'^{t) > 0, 
and 

Step 3: The nth car moves Vn{t) sites ahead. 

The number of empty sites in front of the nth car at time t + 1 can be written as 



For the new model with maximum car velocity I'max — M and a stochastic delay probability 
/, the velocity of the nth car at time step t as a function of the inter-car spacing C„(t) can 
be written as 



Step 1: v'^it) 



mm{Cn{t),M) 



Cnit + 1) = Cnit) + Vn+lit) - Vnit) 



(1) 



Vn{t)^FMU,Cn{t)) 



(2) 



where 



M 



iiC> M 



\ 



C-1 



with probability / 



\ if < C < M 



(3) 



C 



with probability 1 — / 







ifC = 
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III. INTER-CAR SPACINGS IN THE STEADY STATES 



From Eqs.(0)-(|]), we can derive the properties of the inter-car spacings in the steady 
states. Given < M+1, it follows that C„(t) -FM(C„(t)) < 1, and from FM(C„+i(t)) < 
M, we obtain C„(t + 1) = - FA/(C„(t)) + FM{Cn+i{t)) < M + 1. Therefore, if an 

inter-car spacing is not larger than M + 1, it will not be larger than M + 1 as the system 
evolves. 

Given C„(t) > M + 1, it follows that FM(C„(t)) = M, and from FM(C„+i(t)) < M we 
have Cn(t + 1) = — FMiCnit)) + FM^Cn+iit)) < Cnif). Therefore, inter-car spacings 

which are larger than or equal to M -t- 1 will never increase, i.e., if > M + 1, then 

< Cn{t). 

It is useful to define the long and short inter-car spacings via their comparison with the 
maximum car speed M. An inter-car spacing is called a long spacing if it is longer than 
M + 1, i.e., Cn{t) > M + 1. An inter-car spacing is called a short spacing if it is shorter 
than M-|- 1, i.e. C„(t) < M + 1. Based on the above definitions, we can define the excessive 
length of a long spacing Lnif) and the deficient length of a short spacing Snif) as 

L„(t) =max(C„(t)-(M + l),0) 

and 

^„(t) =max((M + l)-C„(t),0). 

It follows that the sum of the excessive lengths of all long spacings L{t) and the sum of 
deficient lengths of all short spacings S{t) are given respectively by 

L{t) = Y,Ln{t), Sit) = J2 Suit) 

n n 

From these definitions, it can be proven readily that 

L{t)-S{t) = Y,[Cn(t)~{M + l)] 

n 

= L-{M + 2) 
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= constant. (4) 

From these properties of the inter-car spacings, we have L„(t + 1) < Ln(t). Hence, 

L{t+l)<L{t). (5) 

From Eqs.(^) and (|^), we have S{t + 1) < S{t). Therefore L and S will never increase as 
the system evolves. If one of the L„ decreases, then L and S will have to decrease. 

Next we look into the question of whether long and short spacings may co-exist in the 
asymptotic steady state. Let Ni{t) be the number of inter-car spacings with length i at time 
t. The probability of finding such a spacing at time t is Pi{t) = Ni{t)/N. Hereafter, Pi{t) is 
denoted by Pi for simplicity, except specified otherwise. Suppose that long and short spacings 
co-exist. Consider a long spacing, if the car ahead moves forward by m — 1 sites, then the 
spacing will decrease by 1. The probability for this to happen is (1 — /)-Pa/-i + /-Pa/. For the 
same reason, the probability for the spacing to be shortened by 2 is (1 — /)Pm-2 + fPhi-i, 
and the probability for the spacing to be shortened by 3 is (1 — f)PM~3 + fPM-2, and so 
on. The probability for the spacing to be shortened by M — 1 is (1 — /)Pi + /P2, and the 
probability for the spacing to be shortened by M is Pq + fPi- On average, a long spacing 
will be shortened by 

M-l 

MPo + 5^ [(1 - /)(M -i) + f{M-z + 1)]P, + /Pm 

1=1 

M-l 

= MPo+Y, - i + /)P, + /Pm (6) 
1=1 

in one time step. The shortened length is positive, unless Pq = Pi = P2 = Pm-i = 

Pm = 0, i.e., S* = 0. Therefore, in the asymptotic steady state of the system, L and S will 
no longer change and at least one of them becomes zero. Hence, it is not possible for long 
and short spacings to co-exist in the asymptotic steady state. 
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IV. ANALYTICAL SOLUTION OF ASYMPTOTIC VELOCITY 

For the low car density case (p < 1/(M + 2)), it is apparent that in the asymptotic 
steady state, L > and S = 0. Hence 

Cn(t)>M + l, Vn. (7) 

In this case, stochastic delay will no longer occur and all the cars will move forward with 
the maximum speed M. The average car speed of traffic flow is 

< V{t oo) >= M. (8) 

For the high car density case (p > 1/(M + 2)), it is apparent that in the asymptotic 
steady state S > and L — 0. Hence 

Cn{t)<M + l, Vn. (9) 

The length of every inter-car spacing cannot be larger than M + 1. Therefore, the average 
speed of traffic flow in the asymptotic steady state is 

M 

<V{t^oo)> = Y^ PSC^ -f) + ii- 1)/] + MPm+1 

i=l 
M 

- J]P,(z-/) + MPm+i (10) 
1=1 

A. vmax = M = 1 

In this case, the high density case refers to p > 1/3, and hence P„ = 0, Vn > 3. It 
implies that only Pq, Pi, P2 are non-zero. To obtain the non- vanishing Pj, we introduce 
Ni^j to describe the number of inter-car spacings with a change in length from i at time t 
to j at time t + 1. The probability of finding an inter-car spacing with length i at time t 
and length j at time t + 1 is 

W,^,{t)^N,^,{t)/N. (11) 
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From Eqs.(|lD-@, we can write down all the non-zero Wi^j as 

Wo^l = Po[{l-f)Pl + P2] 
Wo^2 = 

W^^o = Pi[{l-f)Po + f{l-f)Pi] 

W,^2 = Pl[f{l-f)Pl + fP2] 
W2^0 = 

W2^l = P2{Po + /Pi) 

For |i - j| > 3, Wi^j = 0. 

When the system approaches its asymptotic steady state, all the Pj cease to chang 
the following detailed balance condition for the steady state holds: 

Wi^m = Wm^i, Vm. 

When m equals 0,1, and 2, three detailed balance equations can be written down as: 





.1 + Wo- 


.2 = VTi- 


.0 + W2- 


->o 




.0 + W^i- 


.2 = Wo- 


.1 + W2- 


^1 


W2- 


.0 + W2- 


.1 = Wo- 


.2 + Wi. 


■^2 



Substituting the expressions of Wi^j into the above three equations, we obtain 

P0P2 = /(I - f)Pl 
Normalization requires that Pq, Pi, and P2 satisfy the equations 

J2Pi = P0 + Pl + P2 = l, 
i 

and 

J2 iP^ = Pi + 2P2 = C = 1/p - 1. 



From Eqs.(0),(^4|), and (0), we obtain a quadratic equation for Pq: 

(2/ - ifP^ + [(2/ - ifiC - 2) + l]Po - /(I -f)iC- 2f = 0, (16) 

with its root given by 

p _ -[(2/ - mC - 2) + 1] + v/(2/ - mc -2)C + 1 

n - 2(2/ -1)2 • ^^^^ 



Hence, the asymptotic average speed of traffic flow is 
< V{t ^ oo) > 



C+27^(v/(2/-l)2(C-2)C + l-l) 



-2(-^%^ ^ 2/-1 ^ I 

For / = 1/2, which is an removable singular point, 

< V{t oo) >= C/2 = (1/p - l)/2. (19) 

Equations (|T^) and (19) give the asymptotic < V(t oo) > as a function of / and p in the 
high density case with M = 1. 

B. ^;max = M = 2 

In this case, p > 1/4, hence = 0, Vn > 4. It implies that only Pq, Pi, Pj and P3 are 
non-zero. From Eqs.([l|)-(|^), we can write down the non-zero Wi^j as: 

W^0^1 = Po[(l-/)Pi + /P2] 
Wo^2 = Po[il-f)P2 + P3] 
Wo^3 = 

W,^o = Pi[il-f)Po + fil~f)Pi] 

W,^, = Pi{/(1 - /)Pi + if + (1 - f)']P2 + (1 - /)P3} 
Wi^3 = Pl[il-f)P2 + fP3] 
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W2^0 = P2[il-f)P0 + f{l-f)Pl] 

W2^l = P2{fPo + [f + (1 - fY]Pl + (1 - f)P2} 

W2^3 = P2[f{l-f)P2 + fP3] 

W^^O = 

Ws^, = Ps{Po + fPi) 

W^^2=P3[{l-f)Pl + fP2] 

Substituting the above expressions into the detailed balance condition in Eq.(|12]), we obtain 
the following set of four equations: 

fPoP2 + P0P3 - /(I - f)P? - /(I - f)PlP2 = (20) 

2/P0P2 + P0P3 - 2/(1 - f)P^ - /(I - f)PiP2 - (1 - f)PiP3 + /(I - f)P2 = (21) 

/P0P2 - P0P3 - /(I - f)Pi + /(I - f)PiP2 - 2(1 - f)PiP3 + 2/(1 - /)P| = (22) 

P0P3 - /(I - f)PiP2 + (1 - f)PiP3 - /(I - f)P2 = (23) 



Noting that only two of them, e.g. Eqs. (pO[) and (p3D, are independent. Combining Eqs.(pO| 
and (|23|) with the normalization conditions 



^P, = Po + Pi + P2 + P3 = 1 (24) 

i 

and 

^ «P, = Pi + 2P2 + 3P3 = C = 1/p - 1, (25) 

i 

we can solve for Pq, Pi, P2, P3 and obtain the asymptotic traffic flow velocity for M = 2. 
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V. DISCUSSION 



In order to compare with the analytic results, we carried out numerical simulations on 
ID chain with 1000 cars. The length of the chain was adjusted so as to give the desired car 
density. Periodic boundary condition was imposed. The first 20000 time steps were excluded 
from the averaging procedure so as to remove the transient behavior. The averages were 
taken over the next 80000 time steps. Figures 1 and 2 show the comparison between results 
obtained from numerical simulations and our mean field theory for the cases of M = 1 and 
M — 2 over the entire range of density p. The curves are the theoretical results while the 
symbols represent results of numerical simulations. The curves from the top down along the 
velocity axis correspond to different values of / ranging from to 1. Excellent agreement 
between simulations and our theory is found. 

Prom the fundamental diagrams of the new model, it is noted that when the car density 
is low enough (p < l/(M + 2)), all the inter-car spacings will not be shorter than M + 1, and 
all the cars will not be delayed, leading to traffic flow in its maximum velocity {V — M). 
This situation is more realistic in that in real traffic, no driver would like to slow down his 
car when it is far away from the car ahead. In the high density case, the stochastic delay in 
our new model represents better safety than that of the PI model, and leads to much higher 
asymptotic average velocity of traffic flow than that in the NS model. 

In summary, we introduced a new model with stochastic delays for cars following the 
trail of the car ahead. Its evolution and fundamental diagram are quite different from the 
NS and PI models, even in the simplest case of M = 1. We studied the evolution of the 
inter-car spacings and obtained its fundamental diagram by an analytical COMP approach. 
The results show exact agreement between numerical simulations and our theory. 

The analysis of the dynamical evolution of our new model may give us a clearer physical 
picture on how the acceleration and stochastic delay rules affect the evolution and the 
corresponding asymptotic steady state. It will also provide us with better ideas on developing 
analytical approaches to other traffic flow CA models such as the NS model for which no 
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exact analytic approach has been estabhshed. 
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Figure Captions 



Figure 1: The fundamental diagram with the maximum car velocity M = 1 and for 
different stochastic delay probabilities /. The solid curves are the theoretical results. The 
points with different symbols represent results of numerical simulations. The curves from 
the top down along the velocity axis correspond to different / ranging from / = to / = 1 
in steps of 0.1. 

Figure 2: The fundamental diagram with the maximum car velocity M = 2 and for 
different stochastic delay probabilities /. The sohd curves are the theoretical results. The 
points with different symbols represent numerical simulations. The curves from the top 
down along the velocity axis correspond to different / ranging from / = to / = 1 in steps 
of 0.1. 
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